Non-equilibrium coherence dynamics of a soft boson lattice 
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Abstract 

We study the non-equilibrium evolution of the phase coherence of a Bose-Einstein condensate 
(BEC) in a one dimensional optical lattice, as the lattice is suddenly quenched from an insulating 
to a superfluid state. We observe slowly damped phase coherence oscillations in the regime of 
large filling factor (~100 bosons per site) at a frequency proportional to the generalized Josephson 
frequency. The truncated Wigner approximation (TWA) predicts the frequency of the observed 
oscillations. 



PACS numbers: 03.75.Kk,03.75.Lm 



Proposals for interferometry with BEC have demonstrated that number squeezed states 
can potentially provide robust sub-shot-noise sensitivity to perturbing interactions 
The appeal of one recent proposal by Dunningham and Burnett (DB) is the current 
availability of number squeezed states in an optical lattice 0, |(| and the simplicity of the 
required experimental sequence. In this proposal, an array of number squeezed states is 
initially prepared in a suitably deep optical lattice. The lattice depth is then rapidly low- 
ered and the system is allowed to evolve in the presence of a perturbing potential energy 
gradient. The lattice intensity is subsequently restored to its initial value. A final interfero- 
metric measurement of array phase coherence is used to characterize the perturbation. This 
sequence is loosely analogous to sub-shot-noise optical interferometry using Fock states as 
inputs to a Mach-Zehnder interferometer [lj]. 

In this work, we examine the first stage of the DB interferometer by studying the dia- 
batic response of highly squeezed number states in an optical lattice to a sudden reduction in 
the lattice intensity. This sequence induces on-site phase variance oscillations as the quan- 
tum state evolves. We observe these oscillations through characterization of the dynamic 
evolution of interference contrast. 

Unlike the optical case, analysis of the lattice interferometric sequence is complicated by 
the presence of a strong non-linearity and multiple interfering modes. For the parameters 
of our experiments (many atoms and many lattice sites), exact solutions of the many-body 
equations of motion are unavailable due to the exponentially large Hilbert space needed 
to describe the system. Furthermore, traditional approximations fail due to their inability 
to describe the initial state [Zj. Thus, analysis of this dynamic evolution is an interesting 
problem in its own right. 

The Bose-Hubbard Hamiltonian accurately describes the atom-lattice system 
Written in terms of single particle creation and annihilation operators a] and dj, 

H = -7 53 a\dj + 2^X1 (3id\aldidi + ^ Kajaj, (1) 

<i,j> i i 

where the subscripts index the lowest vibrational mode for each lattice site. Here 7 is the 
tunneling rate between adjacent lattice sites, gfc is the mean field energy due to repulsive 
interactions between two atoms (g = Aiih 2 a/m, a is the s-wave scattering length, and m is 
the atomic mass), Ni the number of atoms per site, and Vi = Qi 2 is the external potential 
due to a harmonic magnetic trap. $ (which depends on the on-site boson occupancy) and 



7 are determined from integrals over single particle wavefunctions The importance of 
quantum fluctuations is determined by the ratio afti/NiPi, where gfii/Nij ~ 1 indicates the 
superfluid to Mott-Insulator (MI) crossover 0, [h], ■ We characterize global array phase 
coherence by the quantity D{t) = J2i^j( a ] a i) /M, where M is defined as the ratio of the 
total number of atoms to the number of atoms in the central lattice site |7(. 

Qualitative understanding of the system dynamics can be obtained by solving Eq. 1 for 
a two lattice site model. Fig. illustrates the evolution of the lattice number distribution 
following a sudden reduction in lattice intensity under conditions similar to those used in 
the experiments described below. Fig. [TJd shows the associated phase distribution (obtained 
through a Fourier transform of the number distribution coefficients). For the two site system, 
the characteristic frequency for these oscillations is determined by the generalized Josephson 
frequency yJANig^ + 47 2 jlH 1^ . At the point of the first phase revival (t ~ 6 ms in Fig. 
lb), the phase variance is sub-Poissonian while the number variance is super- Poissonian. In 
principle, this enables interferometric measurement of phase shifts below the atom shot-noise 
limit p. 

In the case of the full lattice system, we employ the TWA jlfii Q] to obtain approximate 
array dynamics. Under this approximation, the quantum mechanical expectation value of 
an observable is replaced by a semiclassical ensemble average. Specializing to the lattice 
system, we consider a set of wavefunctions ipj which evolve according to the semiclassical, 
discrete Gross-Pitaevskii equation of motion (GPE) HQ: 

= + i> j+ x) + (V, + gf3^ | 2 )^, (2) 

The initial condition ipj(0) = ^Nje l ^ j , where the phases <pj are sampled from a uniformly 
distributed ensemble of values between and 2n, and Nj are determined from the GPE 
groundstate solution. 

The TWA approach is accurate for suitably short times and can be applied to this problem 
since the time evolution occurs in the semiclassical superfluid regime 16j . It has been used 
to analyze the breakdown of adiabaticity for the lattice squeezing experiments described in 



Ref . HI and also to study damping of dipolar motion in an optical lattice [1$, . Density 
matrix renormalization group (DMRG) techniques have recently been used to study non- 
equilibrium dynamics of boson lattices with low filling factor |2(|, but these methods are 



difficult to extend to large filling factors 



2l|. 





FIG. 1: 

a) Evolution of the number distribution of a two site system following a 200 ms lattice intensity 
ramp to 65 Er and a subsequent diabatic lowering of the tunnel barrier to 20 En, for N = 100 
atoms per lattice site. False color exponential scaling denotes probability of number distribution 
in one site, b) Associated phase distribution. 

The apparatus used in this work is described in Refs. on . Rb atoms confined in a 
TOP trap are evaporatively cooled to produce nearly pure condensates in the F = 2, mp = 2 
state which contain ~ 3 x 10 3 atoms at T/T c < 0.3. After adiabatically relaxing the confining 
harmonic potential, the BEC is loaded into a vertically oriented one dimensional optical 
lattice. The lattice is produced by a retro-reflected A = 840 nm laser beam focused to a 
50 yum 1/e intensity radius at the point of overlap with the BEC. This light field provides 
strong transverse confinement in addition to periodic confinement along the propagation 
axis. For example, at a well depth of U — G3Er (Er/Ti = hk 2 /2m ~ 2ttx 3.23 kHz, with 
k = 2tt/X), the transverse oscillation frequency is 150 Hz, significantly larger than the 11 
Hz radial frequency of the magnetic trap [23] . 

We load the atoms into the optical lattice by linearly increasing the lattice intensity to 
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U Q in 200 ms (see Fig. 2a). At U = 63 E R , we infer <y/h = 2tt x 0.019 Hz, N = 90 and 
<?/3o/iVo7 = 2 for the central lattice site (i — 0). We then rapidly lower the lattice intensity 
(in a 2 ms ramp) to a final level Cfy with a corresponding groundstate in the superfluid 
regime. In our experiments, Uf ranges from 12 to 32 Er, covering the parameter intervals 1 
Hz < g(3 /2TTh < 2 Hz and 42 Hz > r y/2irh > 1 Hz. The time scale for lowering the potential 
is chosen to be fast compared to the characteristic time scale for adiabatic evolution of the 
many-body ground state but slow with respect to the oscillation frequency of the individual 
wells. 

We stroboscopically follow the evolution of the array phase coherence after lowering the 
lattice intensity to Uf by holding the atoms in the lattice at Uf before releasing them and 
observing the interference of their de Broglie waves. Fig. |2b displays a series of absorption 
images of the atomic density profile as the quantum state evolves. We define a contrast 
parameter ( as the ratio of the width of a single peak to the separation between the peaks; 
large ( indicates loss of interference contrast. Fig. |2h shows the oscillatory response of ( 
as a function of hold time. Remarkably, the system evolves from an initial state with poor 
coherence to one with sharp coherence and then returns back, possibly indicating a transition 
from an insulator to a superfluid, then back to an insulator. 

To make contact with theory, we study the dependence of the phase variance oscillation 
period T a on the final lattice depth Uf (see Fig. Et). T a is determined through a non-linear 
least squares fits of the functional form exp(— t/r) cos(27ct/T a ) to the oscillation data (r is 
the characteristic damping time and t is the hold time). For 12 < Uf < 32 Er we observe 
a monotonic increase in T a as Uf is increased. The observed period shows a much weaker 
dependence on Uf than, for example, the characteristic period associated with the GPE 



predicted quadrupole breathing mode [24] (solid line in Fig. 3a). 

On the other hand, TWA predicts a weak dependence of T a on Uf, in agreement with 
observations. TWA periods are inferred from the time required for D(t) to reach its maxi- 
mum value (see Fig. Eb) but are scaled by an overall factor of 2.1 in Fig. (dashed line). 
This scaling factor may have its origin in the breakdown of TWA at longer times. In Fig. Efc 
we represent the observed frequencies in terms of effective phonon wave- vectors q e g through 
the (translationally invariant lattice) phonon dispersion relation [25], 



%u q = ,/4 7 sin 2 (^)[2iV^ + 4 7 sin 2 (^)]. (1) 
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FIG. 2: 

a)Experimental sequence for the lattice intensity ramp. The lattice is ramped to a peak value Uo 
near the insulating regime and then rapidly decreased to Uj. b) Absorption images of the atom 
density profile for indicated hold times. c)Phase variance oscillation with Uq = 63 En, Uj = 16.6 En 
and n/h = 2vr x 0.67 Hz. 

Here a = A/2 is the lattice period. We find the value q = q e s such that ui q is equal to the 
observed frequency 27r/T (7 . For sufficiently large Uf the observed frequencies correspond to 
short wavelength excitations, whose frequency is proportional to the generalized Josephson 
frequency. For our parameters, this frequency is much faster than the tunnelling frequency 
which has been observed to determine the time scale for the onset of coherence in lattice 
systems with low filling factor [6] . It is interesting to note that an exponential fit to the data 
extrapolates to 4n/q e fi\ = 1 site at 55 Er, a lattice depth close to the MI cross-over. Fig. Eli 
shows the measured damping coefficient r. As qualitatively expected, faster damping rates 



6 




FIG. 3: 

a) Phase contrast oscillation period vs. lattice intensity. Dashed line depicts results of TWA scaled 
by a factor of 2.1. For comparison, solid line in the upper left hand corner denotes the quadrupole 
mode excitation spectrum calculated using GPE. b) TWA simulation for D(t) with Uf = 23 Eft. 
c) 4-7r/g c ffA lattice sites vs. lattice intensity. Solid line is exponential fit to data, d) Damping 
coefficient r vs. lattice intensity. 



are observed with increased U 



7, 



27| . Future studies of the damping rate as Uf is tuned 



near the superfTuid-MI phase transition may provide interesting insight into the quantum 



critical regime 



26|, |27| . 



Due to non-adiabaticity in the state preparation sequence, we also investigate the possible 
effect of remnant phase coherence in the initial state on the observed dynamics. We do this 
by studying the dependence of the oscillations on Uq, keeping Uf fixed. If the dynamics are 
driven largely by quantum fluctuations, we expect dynamics to be independent of U for 
large Uq (where the associated phase variance approaches 2tt). We find L0 a j2m = 37.6 ± 1.3, 
39.7 ± 1.6 and 38.91 ± 1.6 Hz, for U = 80, 63 and 50 E R respectively with U f = 17 E R 
(Fig. EK). These observations are consistent with the TWA analysis, where D(t) is shown in 
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Fig. for Uq = 80 Er and 50 Er. For comparison, we show numerical solutions to GPE, 
which we obtain by integrating Eqs. 2 over the entire experimental sequence including the 
state preparation, taking (pj = as the initial condition. In contrast to our observations and 
the TWA, these solutions show a strong dependence on the value of Uq (Fig. 0b). Finally, 
we independently assess the role of non-adiabaticity through investigation of exact solutions 
to the two well model for conditions of our experiments. In all cases, we find the number 
fluctuations of the initial state to be at the 1 atom per lattice site level and observe no 
dependence of the observed frequency on residual number fluctuations. 
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FIG. 4: 

a) Phase variance oscillations for Uq = 80 (black), 63 Er (dark grey) and 50 Er (light grey) with 
Uf = 17 Er. b) Numerical model of the lattice sequence using the TWA for Uq = 80 Er (black) 
50 Er (grey) c) Numerical model simulation of D(t) for lattice sequence in Fig. using GPE for 
U = 80 E R (black) and 50 Er (grey). 

In summary, we have shown that the TWA can accurately model the non-equilibrium 
dynamics for the soft boson lattice system. In future work we will seek to study these 
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dynamics in the quantum critical regime as well as to exploit them to demonstrate sub- 
shot-noise interferometry. 
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